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ABSTRACT 



Context. Carbon monoxide is an excellent tracer of the physical conditions of gas in molecular outflows from young stars. 

Aims. To understand the outflow mechanism we need to investigate the origin of the molecular emission and the structure and 

interaction of the outflowing molecular gas. Deriving the physical parameters of the gas will help us to trace and understand the 

various gas components in the flow. 

Methods. We observed CO (12-1 1) line emission at various positions along the LI 157 bipolar outflow with GREAT aboard SOFIA. 

Results. Comparing these new data with CO (2-1), we find basically constant line ratios along the outflow and even at the position 

of the source. These line ratios lead us to estimates of 10^ to 10* cm"^^ for the gas density and 60 to 100 K for the gas temperature of 

the outflowing gas. 

Conclusions. The constrained density and temperature values indicate that we are mostly tracing a low-velocity gas component 

everywhere along the outflow, which is intermediate between the already known cold gas component, which gets entrained into the 

flow, and the hot gas, which gets shocked in the outflow. 
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1. Introduction 

Carbon monoxide is widely considered to be the best probe of 
excitation conditions in the interstellar and protostellar molecu- 
lar gas, because of its high abundance and simple chemistry. In 
the study of outflows from young protostars, in particular, the 
CO molecule plays a central role. Low-7 rotational CO lines 
are routinely used to study the kinematics and energetics of the 
molecular gas swept-out by the protostellar jets and winds. At 
the same time, high-7 rotational CO lines, observable in the far- 
IR wavelength range, are ideal tracers of the gas that is currently 
being shocked in the interaction of the jets with the ambient 
medium, which significantly contribu tes to the overall shocked 
gas cooling (e.g.' Giannini et al]|2001l) . 

In this contribution we report the first observation of a 
higher-7 line of CO, namely the CO (12-11) line, along 
various shocked spots in the LI 157 bipolar outflow that 
were obtained with the GREAT heterodyne spectrometeiQ 
(iHevminck et al.l2012h aboard the Stratospheric Observatory for 
Infrared Astro nomy, SOFIA, dur ing three flights in its Early 
Science phase (lYoung et al.ll2012h . 

The CO (12-11) transition lies at 365 K above the ground 
state and it is therefore a very good tracer of warm molecular gas. 
As such, it can be used to bridge the gap between the low-lying 



' GREAT is a development by the MPI fiir Radioastronomie and 
the KOSMA / Universitat zu Koln, in cooperation with the MPI fiir 
Sonnensystemforschung and the DLR Institut fiir Planetenforschung. 



CO lines, which are related to the swept-out outflow material, 
and the higher-/ lines excited in the shocked hot gas. 

The LI 157 outflow is an ideal target for this pilot study with 
SOFIA because it is one of the most active and extensively stud- 
ied outflows from low-mass young protostars. Observations in 
a wide range of wavelengths indicate the presence of strong 
chemical and physical gradients within the flow. Chemically, 
the outflow is extremely rich, with more than 20 molecular 
species detected, sh owing significant mor phological variations 
along the flow (e.g. i Bachiller et al.l 120011) . This chemical dif- 
ferentiation is associated with the presence of shock spots that 
show a range of excitation conditions that are probed, e.g., 
by multiwavelength H2 observations (Davis & EislofFel 1995]; 
iCaratti o Garatti et al.ll^)06HNisini et al.ll2010ai) . These shocked 
regions show up prominently in H2O emission as well, indi- 
cating that a warm gas-phase chemistry is actively taking place 
(iNisini et al.ll2010bh . 

The outflow has been fully mapped in CO (2-1) by 
IBachiller et al.l (l200lh . while sub-mm CO observations have 
so far only been obtained in the blueshifted southern lobe. In 
particular, Hirano & Taniguchi (2001) have obtained spectra of 
CO (6-5), CO (4-3), and'CO (3-2) in three different shock po- 
sitions, suggesting that the bulk of the emission comes from gas 
atn~ lO'^c m"^ and T ~ 50- 15 K. HERSCHEL observa tions 
obtained bv 'Codel la et all (|2010') and iLefloch et al.l (l2010h sug- 
gested that the CO (5-4) e mission in the b r ightes t blueshifted 
shock spot (called B 1 after IBachiller et al.l (1200 ll) ) can be di- 
vided into different velocity components with different excita- 
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Fig. 1. Positions observed with SOFIA in the LI 157 bipolar out- 
flow, which is emanating from the embedded Class source 
LI 157-mm. The background colour map shows an image in the 
H2 0-0 S(l) Une at n/um (from Nisini et al. 2010b) with su- 
perposed contours o f the integrated CO (2-1) emission (from 
iBachiller et alj I2OOII) . The targeted positions of the presented 
CO (12-11) observations with GREAT are labelled and indi- 
cated with circles, with diameters equal to the actual beam size 
of 21'.'3 of the GREAT beam. 



tion conditions: the low-velocity gas at n ~ 3x10'' cm"^, T ~ 
1 00 K, and the high-velocity gas at n ~ 1 0"* cm"^ and T > 400 K. 
CO transitions at higher J (from J = 15 up to J = 20), which 
probe the hottest gas component, have been detected by ISO at 
low spatial resolution and suggested the prese nce of CO gas 
with temperatures in the range 300 - lOOOK (IGiannini et al.l 
I2OOI !). Very recently, high spatial resolution, non-velocity re- 
solved PACS observations of CO lines from J = 14 to 22 have 
been reported for the LI 157-Bl position, which allows one to lo- 
calize the hot gas in a compact shock of about 9" in size and t o 
better constrain its physical properties (iBenedettini et al.ll2012l) . 

2. Observations and data reduction 

Data of the CO (12-1 1) Hne were obtained with GREAT at var- 
ious positions in LI 157. During SOFIA flights starting out of 
Palmdale, CA, on 14 July 2011 positions on the central source 
and towards the southern outflow lobe were obtained, while on 
the 28 and 30 September 2011 positions towards the northern 
outflow lobe were observed. Fig.[T]shows the six targetted posi- 
tions on a H2 0-0 S(l) and CO (2-1) m ap. They were selecte d 
because they are peaks in the SiO map of IBachiller et al.l(l2001h . 
GREAT was in its Lla/L2 configuration in upper side band. The 
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Fig. 2. Measurements of the CO (12-11) line in the six positions 
in LI 157 marked in Fig.l. Data have been rebinned to velocity 
bins of 3.3 km s"' for these diagrams. The vertical line marks the 
systemic velcocity of visr - 2.7 km s"' . 



fast-Fourier transform spectrometers were used as backends for 
data recording, and data from the XFFTS backend are shown 
here. 

The chop amplitude was 80" to the east and west of the 
source position. The observations were calibrat ed with the chop- 
per wheel method using cold and warm loads (iHevminck et al.L 
2012). Bandpass-averaged zenithal opacities varied from t = 
0.11 to 0.30, with averages of 0.21, 0.14, and 0.17 on the flights 
of 14 July, and 28 and 30, September, respectively. 

The CLASS packagq3 was used for the standard data reduc- 
tion, which included flagging of bad channels and atmospheric 
features, and a third order baseline subtraction. Data in the single 
positions were then averaged using a sigma weighting scheme, 
and calibrated for a forward efficiency of 0.95 and a main beam 
efficiency of 0.54. 

3. Results 

The CO (12-11) line was observed towards the source position 
LI 157-mm and five positions in the low-mass LI 157 outflow 
(Fig.[T]i. The measured line fluxes at the various positions, inte- 
grated from visr - -20kms ' to -1-28 km s' as in IBachiller et al.l 
(I2OOI '), are given in Table[T] CO (12-11) was detected in the 
inner four positions at a signal-to-noise ratio between 3.7 and 
9.5. It was not detected in the northern-most position 1 and only 
marginally, at a signal-to-noise ratio of 2.3, at the southern-most 
position 6 (Fig.|2]l. Although our data do not have at high signal- 
to-noise ratio, we find that at position s 2 and 3 our resolved line 
profiles are resembling those of Bachil ler et all (1200 1 ) in CO (2- 
1), with a peak near the systemic velocity at about 5kms"'. A 
second peak around 10 - 20km s', broader but weaker in the 



see |http://www.iram.fr/lRAMFR/GILDASi 
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Table 1. Measured positions and fluxes. 
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Bachiller et al. data, is seen in our data in these positions as 
well. A similar, but nan^ower peak is also present in position 
5, where a s econdary corn ponent is not seen in the low-J CO 
spectra (e.g. iLefloch et al.l l2010). e obtained the clearest detec- 
tion towards position 5. In this position, which corresponds to 
the bright shocked spot L1157-B1, the emission shows a peak 
close to the systemic velocity of 2.7 km s"', and exhibits a blue 
wing that we can trace to about -20kms"'. 

In their analys i s of lo w- 7 rotational CO emission in LI 157- 
B 1 , iLefloch et al.l (1201 0^ detected two velocity components in 
the lines: a fainter high velocity component (HVC) blueward of 
visr - -7.5 km s"' , and a brighter low velocity component (LVC) 
redward of visr - -7.5 km s"'. Our new observations in CO (12- 
11) are therefore mainly tracing the LVC. 



4. Discussion 

It is interesting to use the new data to estimate the physical con- 
ditions of the emitting gas along the LI 157 outflow, to infer from 
them the origin of the emission, and to see if these conditions are 
changing along the flow. 

To infer the physical conditions giving rise to the CO (12- 
11) emission obs erved with GREAT, we compared it with the 
observations of .Bachiller etal.l (l200ll) . who mapped the entire 
LI 157 outflow in CO (2-1). Their CO (2-1) map was obtained 
at the IRAM 30-m antenna with a beamsize of 11". For the pur- 
pose of our comparison we convolved this map to a beamsize 
of 21'.'3 - con-esponding to that of GREAT for the CO (12-1 1) 
line - and measured the emission in the positions also observed 
with GREAT, integrate d over a velocity ra nge from -20kms ' 
to -1-28 km s', as in Bachiller et al.l (1200 ID . These values are 
reported in Table[T] col. 7. We integrated the emission of the 
CO (12-11) line over the same velocity range (Table[T] col. 6, 
with errors in brackets). The corresponding line ratio of CO (12- 
ll)/CO (2-1) is also reported in Table[T] Its associated error is 
dominated by the rms noise of our CO (12-1 1) observations. 

The first view shows that the CO (12-1 l)/CO (2-1) line ra- 
tio is faMy constant from the source to the innermost northern 
knot LI 157-RO and along the knots in the southern outflow lobe, 
with an average value of about 0.06 + 0.01 . The ratio may be de- 
creasing for the two northern-most knots, but this difference is 
not very significant, given the associated errors, and because for 
the northern knot LI 157-R we can only give a 3-cr upper limit, 
because of the non-detection of CO (12-11). 

It is interesting to note that the CO line ratio of 0.07 observed 
for LI 157-mm, that is "on-source", is not significantly different 
from that of the outflow positions. This suggests that the on- 



source CO emission is also dominated by shock emission, and 
not by emission from the envelope of the embedded source. We 
will investigate this in more detail below. 

To derive information about the physical conditions of 
the CO emitting gas in the observed positions along the 
LI 157 outflow, we used the radiative transfer model RAD EX 
jvan der Tak et al.l 12007 ) to compute the expected CO (12- 
ll)/CO (2-1) ratio to confront it with our observations. We 
adopted a plane-parallel large velocity g radient (LVG) app roxi- 
mation and the coUisional coefficients of lYang et alj (|2010|) . and 
- given the similar flux ratios everywhere - are assuming that 
the emission regions of both lines must be very similar. 

With a single line ratio at hand, it is obvious that the range of 
free parameters has to be restricted with additional information 
to break the otherwise existing strong degeneracies between gas 
density and temperature. To this end, we are considering a col- 
umn density of the CO gas that is compatible with the SPITZER 
observations of pure rotational lines of molecular hydrogen of 
the same regions as we observed here. Nisini et al. (2010b) es- 
timated that the total column density of H2 associated with the 
shocked spots is of the order of 10^° cm"^, with variations of 
only a factor of a few in the various regions. Assuming a CO 
abundance of lO"'*, this translates into a column density of about 
10'^ cm"^ for the CO gas. We also assumed a value for the ve- 
locity dispersion Ay of about 10km s"'. This is what we see in 
our spectrum at position 5 as the width of the line thanks to the 
high spectral resolution of GREAT. 

In Fig. [3] we show (in black) the predicted CO (12- 
1 l)/CO (2-1) line ratio as a function of the kinetic temperature 
and for total H2 densities of lO'*, 10^, and 10^ cm"^. Because the 
line ratios derived from our CO (12-1 1) observations and those 
from the map of Bachiller et al. (2001 ), that are listed in Table[T] 
are so similar for the various positions along the flow, they are 
shown as a black stripe here for clarity. 

Assuming an H2 density of lO** to 10* cm"-', we can con- 
strain the kinetic temperature of the observed CO gas to about 
60 to 200 K. In then' analysis of the single, brightest position 
LI 157-Bl . ILefloch et al.l(l2010l) estimated a gas density of n(H2) 
~ 3 X 10^ cm"^ and a temperature of T ~ lOOK for the gas in 
the low-velocity component. Comparing the estimate from our 
newly observed line ratio with theirs, we find both results to be 
consistent. This confirms that we are mostly tracing the LVC 
component along the whole flow, that is intermediate between 
the cold gas, which gets entrained into the flow, and the hot, 
shocked gas. 

We already mentioned the surprising fact that the CO (12- 
ll)/CO (2-1) line ratio at the position of the central source 
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Fig. 3. Ratios of the integrated fluxes in CO (12-ll)/CO (2- 
1) versus gas temperature for a CO gas column density of 



10 cm at the various observed positions in the LI 157 out- 
flow (in black). For the observed position of the central source 
L1157-mm the ratios for CO (12-ll)/CO (10-9) are given as 
well (in red). 



LI 157-mm is similar to the ratios we are finding along the flow. 
An additional constraint for the CO gas at this central position 
can be provided by an observation of the CO (10-9) line ob- 
tained with the HIFI spectrometer aboard HERSCHEL (Yildiz 
et al., in prep.). This observation was obtained with a beamsize 
of 20", which is comparable to the beamsize of GREAT for the 
CO (12-1 1) line. It gives a flux integrated over the same velocity 
range as our observations of 8.4 + 0.4Kkms"'. This then im- 
plies a CO (12-ll)/CO (10-9) line ratio of 0.4 ± 0.1. In Fig.[3] 
we show the model calculations with RADEX for the line ratio 
of CO (12-1 l)/CO (10-9) (in red) together with that of CO (12- 
n)/CO (2-1) (in black). The assumed column density is again 
10'* cm"^, and gas densities of lO'*, 10^, and lO^'cm"-' are plot- 
ted as lines. The observed line ratios and their associated errors 
are shown as black solid-lined and red dashed bands. 

Fig. [3] shows that a low-density solution would not be able 
to simultaneously account for the two line ratios: the CO (12- 
ll)/CO (2-1) ratio would imply a CO gas temperature of 200 
to 240 K, while the CO (12-1 l)/CO (10-9) would require gas in 
excess of 350 K. Assuming that both line ratios are tracing the 
same bulk mass of gas, densities of between 10^ and 10* cm""* 
and gas temperatures of about 60 to lOOK are required. 

This analysis with the additional use of the CO (12- 
ll)/CO (10-9) line ratio thus confirms and strenghtens the gas 
parameters that we found from the single line ratio along most 
of the LI 157 outflow, and the estimates of Lefloch et al. (201(|) 
for position LI 157-B 1 in the flow. It supports our earlier finding 
that the CO emission in the on-source position is also dominated 



by shocked gas, similar to what we see along the outflow, and 
not by cold emission from an envelope. 

iHirano & Taniguchi (2001) analysed spectra of the CO (6- 
5), CO (4-3), and CO (3-2) lines in three different shock posi- 
tions in LI 157. They found that the emission mostly comes from 
gas at about n ~ lO"* cm"^ and T -50- 150 K. From our SOFIA 
observations we derive a similar temperature range. For the cen- 
tral position, however, a higher density of the order of 10^ to 
10* cm~^ is indicated. 



iBenedettini et al.l ( 1201 2h have analysed the velocity- 
integrated CO emission from J = 14-13 to J = 22-21 observed 
in position Bl (our position 5 ) with HERSCHEL. They derived 
that these lines are compatible with gas at 200 < T < 800 K and 
n > 10^ cm"^. Clearly, CO emission at higher J is dominated by 
a gas warmer than what we are inferring here, testifying that the 
emission of the series of lines of the CO ladder is able to trace 
the temperature stratification that is present across the shocked 
gas in the outflow. 

5. Conclusions 

We presented observations of the CO (12-11) line at various po- 
sitions along the low-mass LI 157 bipolar outflow obtained with 
GREAT aboard SOFIA. 

Comparing these new CO (12-11) data with CO (2-1), we 
found that the line ratios are nearly constant along the flow, even 
at the position of the outflow source. An estimate of the physical 
parameters of the emitting gas shows that we are most likely 
tracing matter with a total gas density of about 10^ to 10* cm"-' 
and at a temperature between 60 to 100 K. 

This indicates that we are mostly tracing a low-velocity gas 
component in the outflow, which is intermediate between the 
cold gas that gets entrained into the outflow and is seen in low- J 
rotational CO lines, and the hot gas in the shocks in the outflow 
that was seen in very high-7 rotational CO Unes with ISO. 
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